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A NOTE ON SPECTRAL PROPERTIES OF THE p-ADIC TREE 


SLAWOMIR KLIMEK, SUMEDHA RATHNAYAKE, AND KAORU SAKAI 


Abstract. We study the spectrum of the operator D*D, where the operator D , introduced 
in [B], is a. forward derivative on the p-adic tree, a weighted rooted tree associated to Z p 
via Michon’s correspondence. We show that the spectrum is closely related to the roots 
of a certain q —hypergeometric function and discuss the analytic continuation of the zeta 
function associated with D*D. 


1. Introduction 

This note builds on our previous paper |6J which described a new spectral triple (A, "H, V) 
for the C*-algebra of continuous functions on the space of p-adic integers Z p . The construc¬ 
tion of this spectral triple utilized a coarse-grained approximation of the space Z p and was 
partially motivated by recent work [3j on spectral triples for more general Cantor sets. Our 
considerations closely resembled standard examples of geometric spectral triples that use the 
usual differentiation for the definition of the operator V. 

The geometric coarse-grained approximation of Z p , which we called the p-adic tree, is a 
weighted rooted tree {V, E} associated to Z p via Michon’s correspondence [3j. The set of 
vertices V of the p-adic tree consist of balls in Z p , with Z p itself being the root of the tree. 
There is an edge between two vertices v and v' if v' C v and v' has the biggest diameter 
smaller than the diameter of v. 

Now consider the Hilbert space H consisting of weighted t 2 functions living on the vertices 
of the p-adic tree: 

H = {/ : V -» C : Y. < <*=}■ (1-1) 

vev 

Here the discrete-valued weight function w : V —> M> 0 is defined by: w(v) = the volume of 
the ball v with respect to the additive Haar measure d p x. For the space of p-adic integers 
Z p , the volume of a ball is equal to its diameter. It is useful to view H as the subspace of 
L 2 (7j p ,d p x) consisting of locally constant functions on Z p . 

Then we introduced an unbounded operator D on H defined on its maximal domain 
{./•f // : /;./>//} bv 
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where deg v is the degree of the vertex v and v' IZ v means that there is an edge between v 
and v 1 . One can think of D as a natural discrete derivative for (complex valued) functions 
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on Z p . This forward tree derivative was then used to construct the Dirac type operator V, 
necessary for the spectral triple. 

It was verified in [6] that the operator D is invertible with compact inverse, implying 
that D*D has compact resolvent. Consequently, the spectrum of D*D is discrete with only 
possible accumulation point at infinity. 

In this paper our main interest is to find the spectrum of the operator D*D. By re¬ 
parametrizing the vertices of the p-adic tree using the set of parameters V = Q v x Z>o, 
where Q p = Q p /Z p is the Priifer p-group, we can decompose H into invariant subspaces: 

where H g = £ 2 (Z> 0 ). This allows the decomposition of the operators D 
and D* into a direct sum of much simpler operators D = ® g eg P D 9’ and D * ~ ®geG P D *g’ 
where D (J is the restriction of D to H g . Identifying the Priifer group with the set of numbers 
{g = : 0 < r < p m ,p \ r}, the operators D g , D* for g = ^ can be written as D g = p m D 0 

and D* = p m D* 0 where D 0 is the operator on £ 2 (Z> 0 ) given by D 0 f n = p n (f n — f n +i)- 
Consequently, D*D = p 2m -DQ.Do and the problem of hireling the spectrum of D*D is 

reduced to finding the spectrum of D^D 0 . 

It will be verified in this paper that the eigenvalues of DqDq are the roots of the g-Bessel 
function i<f > 1 q, A) with q = p~ 2 . In [lj the authors give analytic bounds for these roots 
and discuss their asymptotic behavior. Therefore, we have a good understanding of the 
spectrum of D*D. In particular, using the results of [I], we were able to obtain several 
results on analytic structure and analytic continuation of the zeta function of D*D. 

Part of the motivation for studying the spectrum of D*D is that it may have some relevance 
for developing the structure of p-adic quantum mechanics. The operator D*D, a natural 
analog of the laplacian, can be taken as an alternative starting point for the theory of p-adic 
Schrodinger operators, see [7j. 

The content of this paper is organized as follows. In section 2 we give a brief introduction 
to p-adic harmonic analysis and then describe the p-adic tree associated to the ring Z p . In 
the next section we analyze the forward derivative D on the p-adic tree and its adjoint D* in 
H. We also describe the re-parametrization of the p-adic tree that leads to a decomposition 
of the operator D*D into a direct sum of simpler operators D*D g . Section 4 discusses the 
calculation of the spectrum of DqD 0 along with an “elliptic regularity” theorem that shows 
that the corresponding eigenfunctions are of a special form. In the last section we discuss 
some spectral properties of the operator D*D and the analytic continuation of the zeta 
function associated to it. 


2. Definitions and Notation 

2.1. Fourier Analysis in Z p . We start this section by briefly recalling some of the basic 
results and notation we introduced in [6] regarding harmonic analysis in the space Z p . For 
further reference and more details on this subject we refer to 0,0 and 0. 

The characters on Q p , the set of all p-adic numbers, are given by maps Xa ■ Q P —> C 
defined by Xa( x ) — e 2m ^ ax ^ where a G Q p , and {ax} in the exponent is the fractional part 
of the p-adic number ax. Two characters Xa{ x ) and Xb( x ) are equal on Z p if and only if 
a — b G Z p . Consequently we see that the dual groups of Q p and Z p (denoted Q p , Z p ) are 
Qp = Qp , Z p = Q p /Z p . The dual group of Z p , called the Priifer group Q p , can also be 
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identified with a group of roots of unity given by 

Z p — {e 2 ™^ : n £ Z> 0 ,p \ k £ Z}. (2.1) 

We let £(Z P ) be the space of locally constant functions (test functions), i.e. the set of 
functions 0 : Z p —> C such that for every x £ Z p there is a neighborhood U x of x on which 
0 is constant. The space of linear functionals on £(Z P ) (distributions on Z p ) equipped with 
the weak*-topology is denoted by £*(Z P ). 

If d p x denotes the Haar measure on (Z p , +) normalized so that f z d p x = 1, then we define 
the Fourier transform of a test function 0 £ £ ( Z p ) as the function 0 on Z p given by 

0(H) = / </>{x)Xa(x)dpX. 

J Zp 

For a locally constant function, only a finite number of Fourier coefficients will be nonzero. 
Thus, the Fourier transform gives an isomorphism between £(Z p ) and £(Z p ), where the latter 
in our case is the space of all those functions on Z p that are zero almost everywhere. The 
inverse Fourier transform is given by 

<f>( x ) = Y 0(Hha(H- 

[a] EZ P 

For a distribution T £ £*{Z p ) the Fourier transform is the function T on Z p defined by 
T ( [a ]) = T(xa( x ))- Once again, the distributional Fourier transform gives an isomorphism 
between £*{Z p ) and £*(Z p ). The inverse Fourier transform of a distribution is given by 

T (Xa(x)) = Y T{[a])xa(x). 

[a] GZp 


2.2. The p-adic tree. We recall the construction of the weighted rooted tree {V, E} asso¬ 
ciated to the Cantor metric space (Z p ,p p ), the space of p-adic integers equipped with the 
usual p-adic metric p p , via Michon’s correspondence (3j. [6j. The symbols V and E above 
are used to denote the set of vertices and the edges of the tree, and we call this tree the 
p-adic tree. 

The vertices of the p-adic tree are the balls in Z p . Since Z p is a totally disconnected 
space, the range of p p is countable and consists of numbers of the form p~ n , n £ Z and 
zero. Therefore, if we let V n be the set of balls of diameter p~ n then the set of vertices V 
has the natural decomposition V = IJ^o The set of edges E has the decomposition 
E = (jr=o E n where an edge e = (v,v') between two vertices v,v' belongs to E n if v £ V n , 
v 1 £ V n+ i and v 1 C v. 

The compactness of the space Z p implies that the number of balls of diameter p~ n (hence 
the number of vertices) for fixed n, and the degree of each vertex are finite. Now we observe 
the following fact: 

Proposition 2.1. Every ball of radius p~ n contains a unique integer k such that 0 < k < p n . 
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A proof of this proposition can be found in PJ. From this observation we see that there is 
a one-to-one correspondence between the set of integers 0 < k < p n and the set V n of balls 
of diameter p~ n . Therefore, the set of vertices has the natural parametrization: 

V = S:= {( n , k):n = 0,1, 2,, 0 <k< p n }. (2.2) 

Two vertices (n, k ) and (n + 1, k ') are connected by an edge if and only if k' — k is divisible 
by p n . Thus, a given vertex (n, k) will be connected (via edges) to exactly p vertices in V n+ \. 
Also, we introduce a weight function w : V —>■ M + by w(v) = volume (v) with respect to the 
Haar measure d p x. If v G V n then w(v) = p~ n . 


3. The Operator D 


3.1. A Forward Derivative on the p-adic tree. Due to the decomposition V = IT=0 V n 
of the set of vertices, any complex valued function / on V can be written as a sequence {f n } of 
complex valued functions on V n . Let £*{V) denote the space of all complex valued functions 
living on the vertices of the p-adic tree. By the discussion in the previous section we can 
identify each V n = 7L v lp n 7L v = Z/p”Z with the finite additive group Z/p”Z. Consequently, 
we can introduce the Fourier transform of a function / G £*(V) to be the discrete Fourier 
transform on each V n given by 

1 pn_1 

Inil) = — ^2 fn(k)e~ 2m $, o < l < P n . 

P k= 0 

Because the characters of Z/p n Z satisfy the orthogonality condition: 


— 27 riks 

e p J 

0 <s<pi 



0 if p j \ k 
pi if pi | k, 


we obtain the following Fourier inversion formula: 

fn(k) = ^ /n(0 e pU ■ 

0 <l<p n 


(3.1) 


(3.2) 


We also remark that the p-adic tree can be thought to be self dual, V = V, due to the 
fact that each Z/p n Z (hence V n ) is self dual. Thus, the Fourier transform on the p-adic 
tree is an isomorphism between the space £*(V) of functions on the vertices of the p-adic 
tree and the space 8*{V) of functions on the vertices of the dual tree. Additionally, via the 
Parseval’s identity, the Fourier transform gives an isomorphism between the Hilbert space 
H = £ 2 (V,w), of (11,11) . and H := £ 2 (V), where the latter Hilbert space has no weight in the 
inner product. 

Notice that the decomposition V = (J^o Ki induces a Hilbert space decomposition H = 
Ur^=o^ 2 (^’P _r1 )- Using the parametrization (12. 2 j) introduced in the previous section, we can 
write the action of the operator D of the formula (11.21) on the components f n of / as: 


(Df) n (k) = p n 



1 

P 


Y /n+l ( k + A 0 ”) 

o <j<p 
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We note here that the choice of the domain for D, as well as all other unbounded operators 
below, is the maximal domain, i.e {/ G H : Df G H}. 

Using the Fourier transform of f n and orthogonality of the characters (13. lj) we write the 
following equivalent formula for D\ 


Df„(k) =p n Y, (7»(0 - /»+i(pO) 

0 <1<P 

Hence, in Fourier transform, the operator D becomes D given by: 

Df n (l)=p n (f n (l)-f n+1 (pl)), (3.3) 

which is an unbounded operator on H. Notice that D and I) are unitarily equivalent, but it 
is easier to work with the latter. 

The adjoint D* of D is given by 


D*g n {k)=p ri 


9n(k) - g n -i(k mod p n *) 

P 


assuming p_i(0) = 0. 

Later we will need the following formula for the adjoint D* of D: 


(3.4) 


D*g n (l ) = 


'p n g n {l ) 


if p\l 


P n (dn{l) ~ X p?hi-\ Q)) otherwise. 

It was verihed in [6] that D and D* are invertible with compact inverses. 


(3.5) 


3.2. Invariant Subspaces of H. The key observation that allows us to find the spectrum 
of D*D is that we can decompose the Hilbert space H into invariant subspaces by means of 
a different parametrization of the p-a.clic tree. 

The original parametrization (12.21) of the set of vertices of p-adic tree was done by using 
the set 

S := {( n , k) : n — 0,1, 2,..., 0 < k < p n }. 

Given a pair (n, k) in S notice that we can write k = rp l with p\r and l G {0,1,... n — 1}, 
by factoring out the highest power of p that divides k. Such a representation of k will be 
uniquely determined by r and l. If we associate n with ^ where m — n — l, 

which is a unique representation of n in terms of r and m, then we have the correspondence 

(n,k)^ (^r,z). 

Conversely, given a pair where 0<r<p m ,p{r and l G {0,1,2,...} we can make 

the unique association ^ (m + l,rp l ). Thus, if 


S' 



: 0 < r < p m , p \ r, l 


0 , 1 , 2 ,... 


then we have the one-to-one correspondence between the sets S and S' given by (n, k ) GA 
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In fact, the set of numbers jg = : 0<r<p m ,p\rj is isomorphic to the Pr lifer group 

Q p defined in (12.ip . Therefore, V = V = Q p x Z> 0 . Consequently we obtain the following 
new decomposition of the Hilbert space H\ 


H = f(S) 2* f(S') = ® f? 2 (Z> 0 ) =: ® H g 

^n&Gp g&Gp 

where H g = £ 2 (Z> 0 ). 

We will now look at the operators D and D* in the new coordinates. Using formula (13.31) 
we compute: 



Equation (j3.5j) yields: 

_/ r \ (p m + l f if Z = 0 

\p m J ^ _ if ^ 1- 1 )) otherwise . 

If we assume that / — 1 j =0 for any r, m then we can rewrite the formula for D* as: 


D*f 




(3.6) 


Notice that, in the new coordinates, the operators D and D* affect only the second coordinate 
l and consequently each H g is an invariant subspace. Thus, by letting D g := D\ H and 
D* ■— D*\ Hg of D g , we have the following decompositions of the operators D and D*: 


S = ©U and f = ® S'. 

g(zG P g&Gp 


(3.7) 


Let D 0 be the operator on f? 2 (Z> 0 ) given by D 0 f(l ) = p l ( f(l) — f{l + 1)). It will be more 
convenient to switch to subscript notation and write: 

(■ Dof)n=p n (f n -fn+l )■ (3.8) 

The adjoint of Do is given by 

(D* Q g) n = p n (g n - ^9n-i \ ■ (3.9) 

From formula (13. 6 P we see that if g — ^ then 

D' g = p m Do and D* — p m D^. 

Consequently, D*D has the decomposition; 
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Thus, the key step in finding the spectrum of D*D is to compute the spectrum of the 
operator DqD 0 on £ 2 (Z> 0 ). We devote the next section to a description of this spectrum. 


4. Spectrum of DqD 0 

The fact that D~ l is compact implies that the operators D*D and DqDq have compact re¬ 
solvent. Consequently, the spectrum of the unbounded operator D^D 0 consists of eigenvalues 
diverging to infinity. 

Using formulas (13. 8 j) and (13.9ft we obtain the following system of equations for DqD 0 . 

(AjA)/)o = Zb - h 

{D*D 0 f) n = p n (0of) n - l -(D Q f) n _^j 

= p 2n ~ 2 [-p 2 f n+ 1 + (1 + p 2 )f n ~ fn-i] for any n > 1. 

We remark at this point that we could equivalently study the spectrum of D 0 D J; however 
the equations for the latter operator are not any simpler than the formulas for DqD 0 . Ob¬ 
viously, with the absence of kernels, the eigenvalue equations for both operators yield the 
same eigenvalues. 

The problem is now to solve the following eigenvalue equations for DqD 0 : 


P 2n 2 [-p 2 fn +1 + (1 +P 2 )fn - fn-1 ] = Afor 71 > 1 


(4.1) 


fo ~ fi — A/o, 

with f n e £ 2 (Z> 0 ). 

The key step in solving the system of equations (I4.1ji is the following result which asserts 
that all eigenvectors of D^Dq take the special exponential sum form f n = Yl'k=o c (k)p~ nk with 
rapidly decaying coeffients c(k). This result is a form of elliptic regularity of the operator 

m. 

Theorem 4.1. Let {/ n (A)} be an eigenvector of DqD 0 with eigenvalue X. Then the following 
statements are true. 

(1) The sequence {f n ( A)}£L 0 belongs to £ 1 (Z> 0 ). 

(2) The eigenvector /„( A) can be uniquely expressed in the form 


fn (A) = ^c(2k)p 


—2 nk 


(4.2) 


k =1 


where the coefficients c(2k) decay exponentially in k. 

(3) The coefficients c(2k) satisfy the equations 


c(2) = 


A 


P 


-2 




k =0 


and, for k >2, 
c(2k) = 


-A 


1 — p 


-2 


k -1 


c(2)p k ( k 1 ^ ) ( p 2 — 1) 


k—2 


(' p 4 — l) 2 (p 6 — l) 2 . . . (p 2fc “ 2 — 1 ) 2 fp 2k — 1) 


(4.3) 
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(4) If the remainder r n (2N) is defined by the formula: 

fn( A) = c{2)p~ 2n + c(A)p~ 4n + c(6)p~ 6n + ... + c(2N - 2)p~l 2N ~ 2 > + r n [2N), 
then {r n (21V)}()E 0 —)■ 0 as N —>■ oo in the I 1 norm. 


(4.4) 


Proof. The main idea of the proof is to rewrite the equations (14.11) in an integral equation 
form and then use it iteratively to produce the solution. To this end we regroup the terms 
in the first equation of system (14.11) above to obtain: 

(fn ~ fn- 1) - P 2 (/„+1 - fn) = Xp 2 ~ 2n f n . 

Using the notation A f n := f n+ \ — f n , we can then rewrite the system of equations (14. ip 
as follows. 

A f n = P~ 2 (A/ n _i - A p 2 ~ 2n f n ) for n > 1 

A/o = fi — fo = —A/o- 

Iteratively, with the help of equations if 1.11) . we obtain the following formula for A f n : 

Af n = —p 2n A(/o + /i + ... + f n ), n > 0. (4.5) 

Equation (14.51) is a one-step linear difference equation, so it has one-parameter family of 
solutions. However, since we are looking for the solution in the Hilbert space we need to 
choose one that vanishes at infinity. This leads to the following formula for f n : 

OO l 

fn = Y2 Xp M S f k - 

l—n k= 0 

Interchanging the summation indices of the above formula we obtain: 

n oo oo oo 

fn = Xp ~ 21 f k + E E Xp~ 2l fk 

k= 0 l=n k=n-\- 1 l=k 

A 


(1 ~P~ 2 ) 


—2n 


P 


E A + E P"“ A 


k =0 


k=n -\-1 


(4.6) 


Thus we can estimate: 


y. \fn\ A 

n =0 


A 


(1 -P~ 2 ) 


OO oo 


Ep _ 2 "Eiai + E E p ^ia 


n=0 


k =0 


n =0 k=n-\-l 

By interchanging the summation indices in the first sum above and evaluating the sum 
over n we obtain: 

oo n oo / _ 2 /~ 

Ep _ 2 ”Eiai = Eiai 

n=0 k =0 k =0 


P~ 

1 — p 


-2 


Using Cauchy-Schwartz inequality and the fact that / G £ 2 (N) we conclude that this sum 
is bounded: 

.—2k 


y i a 

k=0 


P 


p- 


, < 

-2 / — 


P 


-2 


P 


-4 


2< OO. 
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Notice that for the second sum we have: 

oo oo oo oo 

£ £ p““iai = Ep“ 2 "Ep“”i/»+iI- 

n =0 k=n +1 n=0 1=1 

Once again using Cauchy-Schwartz inequality we see that the second sum is finite: 


£p- 2 ”£p- 2 UA< 


p 


-4 


n=0 


1=1 


1 — p -4 VI — p 


n -2 


2 < OO. 


This verifies that {f n } E ^(N). 

To prove the second part of Theorem 14.11 we observe that equation (14.6[) gives; 

A p~ 2n 


fn = 


(1 - p ~ 2 ) 


^2fk- Y, f k +J2p 2l f n + i 


_k =0 k=n -\-1 1=1 

Rearranging the terms on the right hand side of the equation to isolate the coefficient of 


p 2n we get: 


fn = 


A 


1 — p 


-2 


£ ft p- 2n - A p- 2n £ 


k =0 


l — p 21 A 


from which we extract the coefficient 

o(2) := 


1 — p 


-2 


1=1 


k =0 


1 — p~ 2 J 


_o I fn+li 


(4.7) 


Notice that c(2) is well deffiied due to part (1). 

Recursively applying this formula for f n on the right hand side of equation (14.711 we obtain: 


fn = c(2)p- 2 " - Ap- 2 " £ 


1=1 


-21 


1 — p 
1 — p~ 2 


c(2)p 


—2n—2l 


A p 


—2n—2l 


£ 

k =1 


1 _ p 2k \ 
1 — p 


) fn+l+k 


Once again we rearrange the terms to extract the coefficient c(4) of p 4n . 

/ OO 

fn = c(2)p~ 2n + I —Ac(2) 


°° 2« _ p-42 


1=1 


—4n 


y2 00 /„-2 1 


P 


An 


£ 

2=1 


This gives: 


c(4) = 


p — p 
1 — p~ 2 

—Ac(2) 


1 — p~ 2 

-42 \ 00 /1 „- 2 fc \ 


P 

1 — p 


£ (•) w 


p 


(1 -p 2 ) (p 4 - 1) 

By repeatedly applying this process we can obtain an expansion of f n in powers of p~ 2n , 
provided the remainder r n (2N) goes to zero as N —> oo. We prove a stronger d 1 estimate on 
r n (2N) below, implying the pointwise convergence needed for the existence of the expansion 

of f n - 

Using induction we readily establish that the coefficients c(2k) of this expansion are in 
general given by the formula: 
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c(2k ) = 


-A 


1 — p 


-2 


k -1 


c(2)p k ( k l \p 2 — 1) 


k—2 


('p 4 — l) 2 (p 6 — l) 2 . . . (p 2fc_2 — l) 2 (p 2fc — 1) 


(4.8) 


Next we estimate the growth of the coefficients c(2k). Simplifying the formula for c(2k) 
we obtain: 


c(2k) = 


c(2)A 


k -1 


2 )( fc !)(1 — yj)(l — ^) 2 (1 — -^) 2 ... (1 — zpr^) 2 (l — 4jr) 


a 1 (i- 


p (-.- 2)(*-.) (1 _ ■ ) 2 (1 _ • ) 2 (1 _ 1 )2 . . . (1 _ 1 ) 2 (1 


P^ 1 / V p^ / \ p u / V p 

rfc /i _-2i\2 \ TT°° 11 2i\2 


~)2 ‘ 


p2k 


Since J} i=:L (l — P 2l ) 2 > n“i(l — P 2 *) 2 an d th e infinite product is a finite nonzero number, 
we obtain the following estimate for the coefficients c(2k): 


\c(2k)\ < 


|A| fc ||/lh 


p(fe—2)(fc—1) — p—2i^2 ’ 

which shows that they decay exponentially. This establishes both the second and the third 
part of the theorem. 

Finally we estimate the remainder term r n (2N). Using induction it is easily established 
that the remainder r n (2N) satisfies the following summation formula: 


r n (2N) 


-A) 


N 


oo oo 


P 


2nN 


££•••£ 

h=l ^2 —1 /at —1 


(p-(2N-2)h _ p -2Nh ) ^p-(2N-4)h _ p -{2N-2)l 2 ^ 


l — p 


-2 


1 — p 


-2 


(1 —P 2lN ) . 

_2 j n-\-i\-\-i 2 -\-.. .-\-i n' 


i— p 

Estimating the i 1 norm we see that: 


|A| 


N 


^ \r n (2N) | < ^2 ^2 ■ ■ ■ ^2 p 2 nN(l _ v - 2 \N 

l 1= ll 2 =l l N =l y y y ' 

...(l-p- 2lN )\f n+h+l2+ ... +lN \. 


,-{2N-2)h _ p -™h\( p -(2JV-4)J a _ p -{2N-2)h^ 


n =1 


Notice that the term E) x5 _ 1 (l — p 2/A, )|/ n+ ; 1+i2+ +i)V | can be estimated as follows: 


(4.9) 


OO oo 

5^(1 - P~ %lN )\fn+h+h+...+l N \ < sup(l -p~ 2lN ) | fn+h+l 2 +...+l N 


< ll/lll 


where the last line can be justified by changing the summation index in the previous line 
appropriately. Moreover, we can explicitly calculate each sum that appears in formula (|4.9jh 
For example: 
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E 

71=1 


P 


-2 Nn 


P 


-2 N 


J2(p-( 2 N-2)h _ p - 2NZi^ = 


h=i 


^(p _(2Ar_4)Z2 - jT (2JV_2) ' 2 ) = 


(2 = 1 


1 — p~ 2N ’ 

(1 — p -2 )p“( 2iV-2 ) 

(l_p-(21V-2))( 1 _ p -21V)’ 

(1 — p~ 2 )p~ ( ' 2N _4 ) 

(1 _ p -(21V-4))( 1 _ p -(2JV-2)) ’ 


and so on. Substituting all these values into the formula (14.9ft we get the following estimate: 


Vl (2N)\<—^— P ^ (1-p 2 )p {2N 2) (1 ~p 2 )p 2 

^ rn[ “ (1 -p~ 2 ) N ' 1 -p~ 2N ' (1 - p-( 2N ~ 2 l)(l - p~ 2N ) ■“ (l-p- 4 )(l-p- 2 ) ' 

Simplifying this expression we obtain: 

(27V) I < _^_<_ WL _ 

n ■ P N(N+v > nL(i-P - 2fe ) 2 ^ (iV+ 1 ) nr=i(i-p- 2fc ) 2 ' 

Since ^^(1 ~~ p~ 2fc ) 2 < oo we see that |r n (2iV)| —> 0 as N —» oo. 

Finally we will prove uniqueness of the expansion of / n (A). Consider the analytic function 
/(z) = From the above estimate of the coefficients c(2k) we see that the 

radius of convergence R of the power series for / is given by: 


= limsup y/|c(2/c)| < |A| limsup 

■ft k —^oo &—>• oo 

Therefore, R = oo and the function /(z) is entire. Therefore, in particular: 

OO 

/(p _2n ) = E c ( 2l '>p' 2nt = aw, 

fc=i 

and so the coefficients c(2k) are uniquely determined by / n (A), because an analytic function 
is completely determined by its values on a convergent sequence of points, [2]. 

□ 


n*i(i-p 


—2z^2 pk-3-\-2/k * 


Remark: The collection of k 2 functions with a power series representation of the form (14.21) 
is fairly restrictive which is clear from the fact that lim„_ 5 . 00 p 2n /(n) = c(2). It can be easily 
shown that the set of k 2 functions with this power series representation is dense in the space 
of all £ 2 functions. 

The difficult part already completed, we can now state our main theorem. 

Theorem 4.2. The spectrum of the operator DqD 0 consists of simple eigenvalues {A n } which 
are the roots of the q-hypergeometric function A ha 1 0 1 q. A), with q = -K. 
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Proof. Substituting f n = Y^k=i c (2 k)p 2nk and formula (14. 8 ft into the initial condition of 
system (14.11) and dividing throughout by c(2) we obtain the following: 


pz 


- (-2k + X _ l \ > k-l 2k-2 

x-i + 22 — - -—— = °- 


fc = 2 n = 2 (i-p 2j ) i- 


,2j—2 


P 


(4.10) 


The inhnite sum on the left hand side of the above equation, call it Si, can be simplified by 
first breaking it up into two terms, extracting some terms and then recombining as follows: 


* = £■ 


(p 2k — i) \ k V 


— 1 — 2 fc —2 


i - 2 R. 2 (i - ?>) (i p2 ,_ 2 


E 


fc—l„ 2 fc—4 


A fc “V 


fc=3 n k Ai-P 2j )[^ 


p 


2 j —2 


£ - 1 ) V 


2 fc—4 \ fc—1 


(1 -P 1 ) (l - Y 


E 

k =3 


p 2k ~ 4 A 


p 2 (p 2k - 1) + (1 - p 2k ) M - 


n- = 2 (i-p 2 j )(i* 


,2i-2 


P 


Using the substitution q — ^ equation (14.10)1 can be written as 


(<?-!) + 


A 


+E 


2—fc \ fc—1 


q 2 ~ k A 


( 1 - 9 ) fen‘d(i-9-i)n‘, 2 (i-9 j - 1 ) 

Notice that at k = 2 the expression 

q.2-k^k-i 

yields the value (1 ^ ) . Thus the above equation is in fact equal to: 


= 0 . 


{q-l) + J2 


2—fc \fc-l 


q 2 ~ k A 


= 0 . 


fc^n^a-rtn^a-^- 1 ) 

Now we rearrange the terms in the inhnite sum in order to compare it with the hyperge¬ 
ometric function i0i q, A). 


(? - !) + HFU 


2—fc\fc-l 


q 2 ~ k A 


(fc- 2 )(fc — 1 ) 


_ , 1) | f- (-') A 'A A ' '(' - q)q^ 

b njacw) 2 


(9 -!) -12 


fc(fe-i) 


(—l) fe A fe (l — 9)9 2 

^ nL (i - Q j y 


i | ” (_i)*A fc g^ 

tt nL(w) 2 


By using the notation 


(a; g) n = (1 - a)(l - aq )... (1 - aq n x ) 
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and the above computation, we can rewrite the eigenvalue equation as 




(-l) fc A k q 


fc(fc-i) 


k =1 


(q\q)t 


o. 


The function i</>i of four variables ao, 6 i, q, z is defined as 


iCl 


a ° ; q 2 , z )=£ 


bi 


( a 0 ;q 2 )r 


^ (g2 ;g 2) n ( 6 i ;? 2) r 


'-i)V®; 


Thus, if A is an eigenvalue, we get: 


m 


°o / i\k \ k fc (fc-!) 

(i; g ,A) = i + X: ( l = 0 , 


fc=l 




showing that the eigenvalues of the operator DqDq are the roots of the above q - hypergeomet¬ 
ric function. Conversely, the above calculation shows that given a root A of A ha 1 (f> 1 q 1 A) 
the formula (14.2p with arbitrary c(2) and other coefficients c(2k) given by (14. 3 p gives, up to 
a constant, the unique eigenvector of DqD 0 corresponding to eigenvalue A. By the analysis 
in jb] the whole spectrum of DqD 0 consists of such eigenvalues. □ 


5. Spectral properties 

5.1. Spectrum of D*D. Computation of the spectrum of D*D is based on decomposition 
(13.1 Op and the analysis of the spectrum of D^D 0 in the previous section. 

Theorem 5.1. Let {A n } be the eigenvalues of the operator DqD 0 and let D*D g be as in 
formula (13 . 7j) . Then, 

(1) The spectrum of D*D g consists of simple eigenvalues {p 2m \ n } i.e., a(D*D g ) = 

U n {P 2mX n}- _ _ 

(2) a(D*D ) = cr(D*D) = U mn {p 2m ^n} ■ Moreover, each eigenvalue of D*D occurs with 
multiplicity p m ( 1 — A). 

Proof. The above results follow directly from the decomposition (I3.10p . Since the number of 
different values of r less than p m that are relatively prime to p is equal to p m — p m ~ l , each 
eigenvalue of D*D in H has multiplicity p m {l — A). □ 

Corollary 5.2. The operator ( D*D )~ 1 is a s-th Schatten class opeartor for all s > 1. 

Proof. From the decomposition (13.101) we see that: 

(D*D)~ S = 0 p~ 2ms ( D*D 0 )~ S , 
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from which we compute the following trace: 

Tr {D*D)~ S = P~ 2 ms Tr(D*D 0 )- s 
■Jm&Gv 

oo 

- £ £ p- 2 "“Ti(D;D„)-. 

m =0 0 <r<p m 
P\r 

Since the number of nonnegative r’s less than p m and relatively prime to p is equal to the 
Euler number of p m , we can compute the sum over m provided that s > \\ 

OO 

Tr(DT>)-' = 

m= ° (5.1) 

= O (t^) 

From JTj] we have that X n < p n , so we can estimate the trace Ti’(.Dq.Do) _s = 
as follows, provided s > 0: 


Tr(BJDo)-” = £(A„r ^ £p“” = T—p- 

n =0 n =0 ^ 

Summing up this information we see that, 

Tr(D D) s < (l - _ pl - 2 ,) 

whenever s > |. Thus for any s > 1 the s-th Schatten norm of ( D*D ) _1 is finite. □ 


5.2. Analytic continuation of the zeta functions. Using formula (15. ip we can express 
the zeta function associated with the operator D*D, denoted Cd{s), in terms of ( Do (s), the 
zeta function associated with the operator DqD 0 : 

o>m = (1 - ^ (rvy^) < 5 ' 2 > 

We now consider the analytic continuation of (d 0 (s). 

Theorem 5.3. ( d 0 ( s ) holomorphic for Res > 0 and can be analytically continued to a 

meromorphic function for Res > —2. 

Proof. To show that (d 0 (s) is holomorphic in the region Res > 0 we estimate: 


l&oMI < £ 


n=l 


^Re s+z Im s 


£ 


\ Re s 
n=l n 


since A n * Ims is unimodular. From [lj we know that the eigenvalues X n of D^Do satisfy the 
following upper and lower bounds: 
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P U 1 


P 


—2 n 


l — p 


—2 n 


< A n < P H . 


Thus we get: 


KdoMi < E 


b" 


Re s ' 


We have an elementary inequality: 

P 2?r 


i — p 


—2 n 


= i__J_>p 2 - 2 


p2n _ p2 _ l ’ 


(5.3) 


which holds since the left-hand side is an increasing function of n, while the right-hand side 
is its value at n = 1. Therefore, we get: 


/„ 2 _i\Res oo 

Eyici. 

which is convergent for Res > 0. Consequently, (d 0 ( s ) is holomorphic in Res > 0. 

We now show that Ct>o( s ) can i 36 analytically continued to Res > —2. Since A„ behaves 
like p n , the analytic continuation of (d 0 (s) will be achieved by a perturbative argument from 
the meromorphic function obtained from the zeta function by replacing A n with p n . First we 
write: 


-ns _ \-s _ -sn lnp _ —sln(A n ) _ 

P A n e e 

Thus, we obtain: 

f—n lnp 


nlnp i 

—e ts dt = s 
dt 


r'—n In p 


e dt. 


' — ln(An) 


- ln(An) 


'‘—nlnp 


\p~ ns — A“ s | < |s| 


e ts \dt = |s| / e mes dt. 


ln(A„) 


J— ln(A„) 

In this integral we can estimate the integrand by its maximum on the interval of integration 
[— ln(A n ), —nlnp] to arrive at the following estimate: 


\p~ ns -K s \< 

Inequality (15,3j) implies that: 

p r 


|s|(nlnp — ln(A n ))e nln P Re ' s if Res < 0 
|s|(nlnp — ln(A n ))e _lll ^ An ^ Res if Res > 0. 


In ( — ) < — In ( 1 — 


P 


—2 n 


—2 n 00 i 

P zn X - 1 


E 


1 — p~ 2n J 1 — p~ 2n j—^ k + 1 \p 2n — 1 


Since 1 — p 2n > ^ for n > 1, we can estimate the above as: 


,n A < 2p -E 


k =0 


k T 1 \p 2 — 1 


= -2 p~ 2n In I 1 


p 2 — 1 
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Consequently, if Res < 0, we have: 

\p~ ns - A" s | < -2|s| In ^1 - P~ 2n P~ nRes 

This lets us estimate the difference of the series as follows: 


n— 1 



< -2 In 



OO 


n =1 


(2+Re s) 


The series p- n ( 2 + Res ) j s convergent for Res > —2 hence, by the Weierstrass M test, 
the series Yl^Li(P~ ns ~ A“ s ) converges uniformly for Res > —2 and hence it is analytic for 
Res > —2. 

Moreover, since 


n =1 


P S 

l — p~ s 


is meromorphic in the complex plane with poles at s = k G Z, we obtain that the zeta 
function A ~ s for meromorphic for Re s > —2 with the above mentioned poles. 


□ 


Corollary 5.4. Cd(s) is meromorphic for Res > —2 with poles at s = and s = 

l (i - ), where k el. 

Proof. The proof of this corollary follows from the theorem above and equation (15.21) . □ 


6. Appendix 


In this section we record some basic properties and identities satisfied by the q - hyperge¬ 
ometric function \fn we encountered in section 4. More on q - hypergeo metric functions can 
be found in |4j. We start with the general definition of these type of functions: 


r+Irs 


CLq ^ Ci \, • • • ^ CL>p 

bi,b 2 , ■ ■ ■ ,b s 


; q, Z 



0 ) q)nifl 1 ) qfn • • • ifln h)n 

(q-, q)n(bi] q)n ■ ■ ■ (b s -,q)n 


(-1 ) n q 


s—r 


Z 


n 


where bj ^ q~ n for any j, n. 

Here we used the notation (a;q) n = (1 — a)(l — aq )... (1 — aq n ~ l ). We remark that 
(a; q)oo = R°0 O (1 — o.q^). When s > r the above series converges for all z while it converges 
for \z\ < 1 when s = r. 

We are interested in the special case where r = 0, s = 1 and a = 0, b = q, which leads to 
the formula: 


0 

q 


q, z = 


E 




-(-i ) n q yp z r ‘ 


The function ifi satisfies the Cauchy’s sum: 


i0i 


a 

b 


; q, b/a 


(b/ ct; q)oo 

(b; q ) 0 o 
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In |T] the authors investigated the roots of the third Jackson g-Bessel function: 

7(3) 1 \ x ( 0 

Ju q) ■= z —7— i 0 i n v +1 ; q, qz 

( 9 ; 9 ) 0 o V 9 

where 0 < q < 1 and z is a complex parameter. It is known that this function has infinitely 
many zeros, each of multiplicity one, all of them real. When v = 0 we that the third Jackson 
g-Bessel function equals the function ifa which we used in this paper. 

We record the following transformation property of 



101 


0 

b 


q, z = 


(z;g)c 

0; q) c 


m 


; q, b 


Starting with this transformation the authors in [dj deduce that if q < (1 — q) 2 then the 
positive roots 04 (g), k — 1, 2, 3,..., of Jq 3 \z; q ), arranged in the increasing order satisfy the 
following: 


q -k/2+a k ( q ) K < g -k/ 2 ? 


where 


oik(q) = 


log 1 


1 -q k 


log q 


In particular, this gives the asymptotic behavior 04 ~ q~ k / 2 as k —* 00 . 
results give an upper and lower bound (15.31) for the roots of the specific 
in this paper. 


Additionally, those 
! 0 ! function needed 
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